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1  Introduction 

The  question  of  which  factors  contribute  to  mixing  in  the  ocean  has  been  the  subject  of 
many  studies.  The  idea  that  organisms  could  play  a  role  in  affecting  their  environment 
through  something  like  mixing  first  appeared  in  Darwin’s  last  book,  called  “The  formation 
of  vegetable  mould  through  the  action  of  earth  worms,  with  observations  on  their  habits.” 
(See  figure  1.)  This  was  suggested  by  his  uncle  and  future  father-in-law  Josiah  Wedgwood 
II,  son  of  the  famous  potter.  From  Darwin’s  book: 

In  the  year  1837,  a  short  paper  was  read  by  me  before  the  Geological  Society 
of  London,  “On  the  Formation  of  Mould,”  in  which  it  was  shown  that  small 
fragments  of  burnt  marl,  cinders  &c  which  had  been  thickly  strewed  over  the 
surface  of  several  meadows,  were  found  after  a  few  years  lying  at  the  depth  of 
some  inches  beneath  the  turf,  but  still  forming  a  layer.  This  apparent  sinking  of 
superficial  bodies  is  due,  as  was  first  suggested  to  me  by  Mr.  Wedgwood  of  Maer 
Hall  in  Staffordshire,  to  the  large  quantity  of  fine  earth  continually  brought  up 
to  the  surface  by  worms  in  the  form  of  castings.  These  castings  are  sooner  or 
later  spread  out  and  cover  up  any  object  left  on  the  surface.  I  was  thus  led  to 
conclude  that  all  the  vegetable  mould  over  the  whole  country  has  passed  many 
times  through,  and  will  again  pass  many  times  through,  the  intestinal  canals  of 
worms,  [emphasis  added]  Hence  the  term  “animal  mould”  would  be  in  some 
more  appropriate  than  that  used  of  “vegetable  mould.” 

The  modern  name  for  this  phenomenon  is  ‘bioturbation.’  At  the  time,  the  idea  that  many 
small  events  could  accumulate  to  create  something  much  bigger  was  quite  controversial  (and 
still  is  to  some  people,  in  the  case  of  evolution). 

It  was  not  until  almost  a  century  later  that  Walter  Munk  [8]  seriously  proposed  marine 
organisms  as  a  possible  factor  for  mixing  in  the  ocean.  However,  Munk  discounted  the 
effect  as  negligible.  The  idea  lay  dormant  for  forty  years,  until  2004  when  Huntley  & 
Zhou  [3]  analyzed  the  swimming  of  100  species  in  a  wide  range  of  sizes,  and  concluded 
that  for  11  representative  species  the  turbulent  energy  production  is  10  5  W  kg"1,  the 
total  of  which  would  add  up  to  an  amount  comparable  to  the  energy  dissipation  by  major 
storms.  Another  study  on  the  magnitude  of  the  mixing  effect  of  the  oceanic  biosphere  was 
performed  by  Dewar  et  al.  [2],  who  estimated  that  the  ocean’s  organisms  take  in  a  total 
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Figure  1:  Darwin’s  last  book,  “The  formation  of  vegetable  mould  through  the  action  of 
earth  worms,  with  observations  on  their  habits.” 


power  input  of  almost  63  TeraW.  They  assumed  that  about  1%  of  this  amount  is  delivered 
as  mechanical  energy,  which  leads  to  a  number  comparable  to  the  input  of  winds  and  tides 
[2].  An  experimental  study  by  Kunze  et  al.  [5]  showed  that  the  level  of  turbulent  activity  in 
an  inlet  was  elevated  by  2  to  3  orders  of  magnitude  during  the  day,  due  to  swimming  krill. 
(Though  it  must  be  mentioned  that  later  studies  were  far  less  successful.) 

In  reaction  to  the  papers  mentioned  above,  Visser  [10]  argues  that  very  small  organisms 
such  as  zooplankton  cannot  produce  sufficient  mixing  to  overturn  a  stratified  medium.  The 
author  points  out  that  the  efficiency  of  mixing  is  proportional  to  the  ratio  of  the  lengthscale 
of  the  organism  (L)  and  a  buoyancy  lengthscale,  defined  as  B  =  ( e/N 3)1^2,  where  e  is  the 
rate  of  turbulent  energy  dissipation  and  N  is  the  buoyancy  frequency.  When  L  <  B,  the 
mixing  efficiency  quickly  drops  by  several  orders  of  magnitude  [10],  rendering  biomixing 
negligible.  Katija  and  Dabiri  [4]  however  suggest  that  it  is  not  the  scale  of  turbulence 
produced  by  the  marine  organisms  that  is  the  relevant  quantity,  but  the  net  displacement 
of  fluid  particles  they  cause.  This  aspect  will  be  discussed  in  these  notes  by  constructing  a 
simple  model  of  stirring  by  swimming  organisms,  based  on  drift  trajectories  due  to  a  moving 
object. 

2  Displacement  of  a  fluid  particle  due  to  a  moving  body 

We  will  first  discuss  the  displacement  of  a  fluid  particle  due  to  a  single  moving  body,  and 
generalize  this  theory  in  the  next  section  to  multiple  ’kicks’  due  to  several  moving  bodies, 
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and  ultimately  to  an  effective  diffusivity.  For  simplicity  we  will  assume  the  moving  body 
to  be  2D  and  symmetric  with  respect  to  its  direction  of  motion.  (Axisymmetric  about  its 
direction  of  motion  in  3D.)  (See  figure  2.)  We  only  consider  potential  flow.  The  body 
moves  along  a  straight  line  with  constant  velocity  U  for  distance  A.  The  distance  between 
the  fluid  particle  and  the  straight  line  of  motion  is  a,  while  the  distance  from  the  start  of 
the  trajectory  to  the  point  of  closest  approach  (with  respect  to  the  initial  position  of  the 
particle)  is  b.  Together  a  and  b  are  called  impact  parameters.  Note  that  a  is  positive  but 
b  can  have  either  sign,  with  a  negative  b  indicating  a  swimmer  that  begins  its  trajectory 
already  having  passed  the  fluid  particle. 


Figure  2:  Displacement  of  a  fluid  particle  by  a  swimming  body. 


Due  to  the  movement  of  the  swimming  body,  the  fluid  particle  will  have  a  net  displace¬ 
ment  A  =  A\(a,b),  known  as  drift.  Although  Maxwell  [7]  was  the  first  to  discuss  the  drift 
to  a  moving  body,  it  is  often  referred  as  ’Darwinian  drift’  [1],  due  to  Charles  G.  Darwin, 
great  grandson  of  the  other  Darwin. 

We  first  need  to  determine  how  much  a  fluid  particle  is  displaced.  We  could  solve  this 
problem  in  a  fixed  frame  centered  around  the  moving  cylinder,  where  the  particle  moves 
around  the  body  (figure  3,  left  panel).  It  is  however  easier  to  work  in  a  frame  which  is 
moving  with  the  fluid  particle  (figure  3,  right  panel). 

In  the  x  direction,  the  displacement  of  the  particle  solely  due  the  flow  U  would  be  b  —  A 
(“free-streaming”).  Because  of  the  influence  of  the  moving  body,  the  fluid  particle  ends  up 
at  X{,  a  distance  Ax  away  from  b  —  A.  Now  we  can  define  a  travel  time  T  it  takes  the  fluid 
particle  to  go  from  b  to  xf 


T  _*_r  * 

U  Jh  u(x,  y) 

where  u(x,  y )  is  the  velocity  felt  by  the  fluid  particle, 
component  of  the  velocity  field  never  vanishes.  Using  Xf 


(1) 


This  assumes  that  the  horizontal 
=  b  —  A  +  Ax,  this  gives 


a  _  rb~x+ Ax  dx 

U~  Jh  v 


c  dx 

M 

dx 


[b  — 
lb—  \+Ax  1^1 
rb-\ 


If  the  particle  is  only  lightly  displaced  and  \b  —  A|  is  large,  then  the  velocity  felt  by  the 
particle  approximates  the  free-stream  velocity  (|u|  ~  U )  in  the  second  integral.  The  second 
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FIXED  FRAME 


COMOVING  FRAME 


Figure  3:  Left:  displacement  of  the  particle  in  a  fixed  frame.  The  particles  starts  in  the 
origin.  Right:  displacement  of  the  particle  in  a  frame  that  moves  with  the  particle.  The 
moving  cylinder  is  in  the  origin. 


integral  can  then  be  approximated  by  A x/U,  which  gives 

A  fb  dx  Ax  fb  dx  A 

U  ~  Jb_x  \u\  U  Jb_x  \u\  U 


From  this  statement  we  can  derive  an  equation  for  the  approximate  horizontal  displacement 
of  the  fluid  particle  Ax: 


Ax  ~  U 


(2) 


The  displacement  in  the  y-direction  is  trivial,  as  the  fluid  particle  will  follow  a  streamline 
of  the  flow  in  the  comoving  frame  (assuming  the  flow  is  steady  in  that  frame).  Far  away 
from  the  moving  cylinder  (A  — >  oo,  b  —  A  -A  — oo)  the  displacement  in  the  y-direction  must 
therefore  be  zero: 


Ay  =  0 


(3) 


Note  that  both  boxed  expressions  are  only  valid  in  the  limit  of  large  path  length  A. 


2.1  Example:  A  cylinder  in  potential  flow 

As  an  example,  consider  a  cylinder  in  potential  flow  with  stream  function 

*(x,y)  =  -Uy(  1-^rV)  (4) 

V  £  +y  / 

with  l  the  radius  of  the  cylinder  and 

ST  ST 

U  =  Sy  ’  V  =  ~dx' 

We  will  choose  U  =  1  and  1  =  1. 
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2.1.1  Far  away  from  the  cylinder 


It  can  be  expected  that  the  effect  of  a  moving  cylinder  on  a  particle  far  away  from  the  body 
is  negligible.  To  test  this  statement,  one  can  construct  a  qualitative  order  of  magnitude  for 
the  displacement  of  such  a  particle.  This  displacement  can  be  represented  by  Ax  and  Ay. 
We  will  assume  the  path  length  is  infinite,  so  the  displacement  Ay  vanishes  identically. 

First  of  all,  let  us  calculate  the  maximum  excursion  in  y,  emax.  This  can  be  calculated 
using  the  stream  function  in  the  fixed  frame  (figure  2).  Here,  the  particle  is  located  initially 
at  (6,  a).  Because  the  flow  around  the  cylinder  is  steady,  the  particle  moves  along  the  stream 
line.  Therefore,  we  can  set  up 

y  a  . 

~y+^T, - 2  =  +  2  ,  U2  -  ~Q’  (5) 

xz  +  y  az  +  bz 

where  a  3>  1  is  assumed  such  that  can  be  neglected.  Let  y  =  a  +  e,  where  a  3>  e.  If 

we  input  this  expression  into  equation  5,  we  obtain 


—a  —  e  + 


a  +  e 


x2  +  (a  +  e)' 


~  —  a  —  e  + 


x2  +  a2 


~  —a . 


Therefore, 

Here,  emax  ~  maxx(e),  thus 


e  ~ 


x2  +  a 2 


0{l/a) 


(6) 

(7) 


where  we  assume  from  now  that  x  <C  a. 

Ax  can  be  calculated  from  equation  2.  First,  we  can  find  u  using  the  given  streamfunc- 
tion  in  the  fixed  frame  and  then  construct  an  approximate  expression  for  |rt|_1  —  1.  The 
velocity  u  is 

dT  x2  —  y2 

11  —  ,  . -  =  -  1  +  - Q  • 

9y  ( x 2  +  y2) 


Therefore,  using  i>i/>  1,  it  follows  that 
1 


-  1  = 


u 


1  - 


1  x2,  —  v2  x2  —  v2 

+  — — *U-i=  x  y 


x*—y* 

(x2-\-y2)2 


( X 2  +  y 


2\2 


(x2  +  y 


2^' 


Let  y  =  a  +  e,  this  gives 


-  1 


u 


x2  —  a2 

t2~\2 


2ae 


(: x 2  +  a2)  (x2  +  a2) 


2^2- 


Therefore, 


r°°  f  i 

J — oo  \  1^1 


Ax  =  /  —  1  )  dx  ~ 


[°°  x2  —  a2  ,  2  ae 

/  - 7T  dx  —  /  - o  dx. 

Loo  ( X 2  +  a2)2  Loo  (x2  +  a2)2 


Here,  the  first  integral  on  the  right  hand  side  is 


/ 


oo  2  _  2  Iff 

dx  =  — r  cos2dde  =  o, 


-oo  (x2  +  a2)^ 


a  J-l 
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where  x  =  atan0  is  used.  The  second  integral  on  the  right  hand  side  is 


2  ae 


L 


-oo  (x2  +  a2) 


2  ae 

2dx  =  ^r 


cos2  6d6  = 

a1 


Therefore, 


Ax  — 


(8) 


Previously,  we  showed  that  e  ~  0(a  x),  which  leads  to  the  conclusion  that  Ax  0(a~ 3). 
Because  a  is  much  larger  than  1 ,  the  displacement  of  the  particle  far  away  from  the  cylinder 
is  negligible,  which  confirms  our  initial  guess.  However,  note  that  the  perturbation  to  the 
velocity  field  due  to  the  moving  object  is  of  order  a-2,  and  the  net  displacement  is  smaller 
by  another  order  of  magnitude.  Thus,  drift  trajectories  far  from  the  moving  body  are  very 
nearly  closed. 


2.1.2  Close  to  the  cylinder 

A  more  interesting  limit  is  the  one  very  close  to  the  cylinder  (a  <C  1).  The  path  of  the  fluid 
particle  can  be  subdivided  into  three  regions  (figure  4).  In  the  first  region  the  particle  more 
or  less  follows  the  x-axis.  In  region  2  the  fluid  trajectory  bends  upwards  and  in  region  3  it 
follows  an  almost  circular  path  to  the  y-axis.  On  the  other  side  of  the  moving  body,  the 
fluid  particle  follows  the  same  path  in  reverse  order. 


Figure  4:  Path  of  a  fluid  particle  close  to  the  cylinder  ( a  -C  1). 


We  want  to  know  the  displacement  Ax  (equation  2).  With  U  =  1  the  integral  for  Ax 
becomes 

T =/ (;  +  1)'iI 


where  we  now  use  the  travel  time  T  instead  of  distance  Ax  to  match  the  units  (T  has  the 
same  numerical  value  as  Ax  for  U  =  1).  We  can  split  the  integral  into  the  three  regions. 
For  the  third  region  one  can  use  the  full  stream  function,  but  for  the  first  two  it  is  easier  to 
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(9a) 

(9b) 

(9c) 


Region  2 

For  region  2  we  first  need  to  specify  the  boundaries  of  the  region.  Let  us  rewrite  the 
approximate  streamfunction  as 

T2  —  -2  (x  -  1)  y  =  -2 XY 

where  X  =  x  —  1  and  Y  =  y.  The  fluid  particle  enters  region  2  at  (Xq.  Lo).  At  that  point 

=  —2XqYq  =  —  (l  —  b~2)  a . 


As  Xq  =  e  and  b  3>  1,  this  is  approximately 

-2elo  =  -a  Y0  =  . 

Because  the  stream  function  is  hyperbolic  and  the  domain  of  region  2  is  of  size  e  by  e,  it 
follows  that  (Xi,Yi)  =  (a/2e,  e).  The  travel  time  in  region  2  is  thus 


T2 


(i  + 1)  dx 


dx 
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Figure  5:  Close-up  of  the  path  of  the  fluid  particle  in  region  3. 


For  the  displacement  of  the  fluid  particle  we  need  the  total  travel  time,  which  is  twice 
the  sum  of  the  travel  time  in  region  1,2,  and  3.  Adding  up  the  leading  orders  of  the  travel 
times  gives 

T  =  T,  +  T2  +  T3  =  (l  log  (?)  -  6-1)  +  (-1  log  (^))  +  (-1  +  1  log2  -  1  loge)  . 

The  e’s  cancel  out.  To  leading  order,  the  total  transit  time  for  half  the  particle  path  is  then 

T  =  —  \  logo  -  1  +  |  log 2  —  b~l  .  (10) 

For  small  a,  the  first  term  dominates  the  travel  time,  i.e.,  the  particle  tends  to  get  stuck 
near  the  stagnation  point  in  region  2  for  a  while.  As  a  — >  0,  we  see  that  T  — >  oo.  In  general, 
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for  a  moving  body  of  any  shape  in  2D  potential  flow  it  turns  out  that  the  coefficient  of  log  a 
is  given  by  the  sum  over  the  linearization  coefficients  for  each  hyperbolic  stagnation  point 
encountered,  as  long  as  a  is  small.  The  displacement  of  a  fluid  particle  close  to  a  moving 
body  in  potential  flow  due  to  this  moving  body  is  thus  well-approximated  by 


Aa  (a,  b) 


—  log  a  0  <  b  <  A 
0  otherwise. 


(11) 


The  cutoff  for  b  <  0  and  6  >  A  is  due  to  the  fact  that  since  the  displacement  is  mostly  due 
to  the  stagnation  points,  the  particle  must  at  least  reach  the  stagnation  point.  (We  take 
the  size  of  the  body  itself  as  negligible.) 

3  Displacement  of  a  fluid  particle  due  to  multiple  kicks 

In  previous  section  an  expression  for  Aa  (a,  b)  was  obtained.  Based  on  this  information  we 
can  derive  an  effective  diffusivity.  Here,  we  have  four  constants  and  two  random  variables: 
U  is  the  mean  velocity  of  the  cylinder,  l  is  the  length  scale,  A  is  the  mean  free  path  of 
the  cylinder,  and  n  the  number  density.  The  position  of  the  particle  is  coordinated  by  two 
random  variables,  a  and  b,  which  were  described  in  the  previous  section.  If  we  pick  a  point, 
(x,y),  in  Cartesian  coordinates,  the  probability  measure  of  the  point  is  simply  represented 
as  (1  /V)dxdy,  where  V  is  the  total  volume.  Instead  of  (x,  y)  space,  we  have  to  use  (a,  6) 
space  for  representing  the  PDF  of  the  particle  distribution.  Considering  that  a  is  positive 
and  neglecting  the  edge  of  the  domain,  we  find 

1  2 

—  dxdy  =  —  dadb .  (12) 

A  similar  relation  can  be  found  in  3D.  In  this  case,  the  probability  measure  can  be  expressed 

aS  12 

—  dxdydz  =  dadb .  (13) 

Now,  assume  that  our  target  particle  is  kicked  by  a  swimmer.  For  N  steps  or  kicks, 

N 

XN  =  x0  +  a  {ak,bk)rk,  (14) 

k=l 

where  ak,  b/~,  and  rp-  are  random  variables.  Their  choices  at  each  step  are  independent.  In 
particular,  (fk'^l)  =  Sp;p.  The  direction  of  a  swimmer  is  chosen  independently.  Without  loss 
of  generality,  we  can  assume  that  xq  =  0.  Since  (rk)  =  0,  it  follows  from  equation  (14)  that 
(xn)  =  0.  To  calculate  the  effective  diffusivity,  the  second  moment  must  be  considered: 

N 

(|*Ar|2)  =  £<A  («fc)  bk))  +  vanishing  cross  terms. 
k= l 

The  cross  terms  vanish  because  {ry.  ■  r{)  =  5ki ■  Each  step  of  the  target  particle  caused  by 
the  swimmers  is  independent  such  that  the  particle  movement  can  be  considered  to  be  a 
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random  walk.  Since  the  variables  are  identically  distributed,  we  get 


v— ^  N  f  r 

(ak,  bk)  hh)  =  N(AX  (a,  b))  =  —  /  /  A2  (<b  b)  2 dadb. 

k= 1 


N  is  introduced  as  the  number  of  time  steps.  Physically,  N  can  be  understood  as  the 
number  of  collisions  between  the  particle  and  the  swimmers.  It  can  be  represented  as  t/T, 
where  t  is  an  elapsed  time  and  T  is  the  mean  free  time,  the  average  time  for  a  swimmer  to 
collide  with  the  particle.  Hence,  t/T  is  the  average  number  of  collisions  during  the  time  t. 
Here,  T  can  be  represented  as  A /U,  where  A  is  the  mean  length  which  a  swimmer  proceeds 
before  changing  its  direction,  so  that  N  can  be  written  as  Ut/X.  Hence, 

(|aV)|2>  =  YV  JJ  (a'h^2dadb' 

where  x(t)  is  used  instead  of  xn-  To  eliminate  V  from  the  equation,  observe  that  the  above 
equation  is  for  one  swimmer,  so  n  =  1/V  is  the  number  density.  Therefore, 

( |m(t) | 2)  =  - -j—  J j  A2x(a,b)  dadb  =  Ant,  (15) 

where  k  is  the  effective  diffusivity.  Ant  comes  from  the  result  of  random  walk  (see  previous 
lecture).  The  above  expression  is  for  a  two-dimensional  system,  but  can  be  easily  generalized 
for  3D.  We  obtain  finally 


(16) 


Recall  our  approximate  form  for  displacement  due  to  a  cylinder,  including  dimensions: 


f  TX  I  Aa  (a>  b)  dadb, 

2D 

H\j  -  \ 

\  TUX  I  tX2x(a,b)  adadb, 

3  D 

(  —l\og(a/l)  0<b<\ 

\  negligible  otherwise  ’ 


the  calculation  of  k  becomes 

Un 

k  ~  — —  A  /  l2  log2  (a /l) da  . 

2A  J  o 

Using  fj  log2  xdx  =  2,  we  obtain 

k  ~  Uni3  .  (17) 

Note  that  this  result  is  completely  independent  of  A  (for  A  large). 

Another  relevant  case  is  a  swimmer  with  a  ‘bubble  wake,’  that  is,  a  region  enclosed 
by  streamlines  that  follows  the  swimmer,  often  called  the  atmosphere  of  the  swimmer.  If 
there  is  a  bubble  wake  behind  a  swimmer,  a  particle  inside  the  bubble  wake  will  follow  the 
swimmer  until  the  swimmer  changes  its  direction.  Therefore,  the  particle  trapped  inside 
the  bubble  moves  by  A.  Our  approximate  expression  of  Aa  (a,  b)  is 


f  A  0  <  b  <  A 
(  0  otherwise. 


(18) 
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Swimmer 

A-dependence 

far/near  field  dominance 

potential  (slip) 

none 

near 

viscous  (squirmer) 

none 

far 

viscous  (no-slip) 

log  A 

near 

trapped  bubble 

A 

near 

Table  1:  Dependence  of  displacement  on  path  length  and  whether  the  displacement  is 
dominated  by  far  or  near  field,  for  various  swimmer  models. 


In  this  case,  our  calculation  of  n  is 

6k  =  ‘—r—  [  \2dadb  =  UnAVbubbic  > 

^  J  bubble  interior 

where  Vbubbie  is  the  total  volume  of  the  bubble.  Therefore, 

k  =  gDnAVbubbie  •  (19) 

Now,  unlike  the  potential  flow  case,  the  effective  diffusivity  depends  explictly  on  path  length 
A.  This  k  can  be  much  larger  than  that  for  untrapped  fluid. 

More  complexity  could  come  in  if  we  consider  viscous  swimmers  with  a  boundary  layer. 
In  this  case,  some  simple  estimates  [9]  suggest  k  ~  log  A.  For  micro-organisms  in  the 
Stokes  flow  approximation,  the  transport  is  larger  and  is  dominated  by  far-field  (a  few 
body  lengths)  hydrodynamics  [6].  The  dependence  on  flow  characteristics  and  boundary 
conditions  is  summarized  in  Table  1,  but  it  should  be  noted  that  the  details  of  this  table 
are  still  the  object  of  active  research. 
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